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A direct sequence code division multi- 
ple access (DS-CDMA) receiver comprises an 
adaptive filter (8) controlled by an adaptive al- 
gorithm (10) for filtering data which has been 
multiplied at (2) by a spreading code, the filter 
having a length equal to the number of chips in 
the code, and a multiuser detector (14) operat- 
ing on the output of the adaptive filter. Prefer- 
ably, either the fast a-posteriori error sequential 
technique (FAEST) algorithm or the stabilised 
FAEST (SFAEST) algorithm is used as the al- 
gorithm (10). 
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ADAPTIVE FILTER FOR A DIRECT SEQUENCE SPREAD SPECTRUM RECEIVER 

Background to the Invention 

The present invention relates to a telecommunications 
receiver employing a new direct sequence code division 
5 multiple access (DS-CDMA) architecture which allows the use 
of fast adaptive algorithms. 

Two adaptive algorithms are commonly in use, the LMS 
and RLS Algorithms 1 ' 2,3 and these are described in Appendix I. 

The least mean square (LMS) algorithm (and the closely 
10 related normalised least mean squares (NLMS) algorithm) is 
a stochastic gradient algorithm which has only one 
parameter, the step size /i. The LMS algorithm is 
computationally simple but its convergence rate is slow and 
highly dependent on the properties of the input signal, more 
15 specifically on the eigenvalue ratio of the autocorrelation 
matrix. When many elements of the input signal are unknown, 
for example the channel in a mobile communications system, 
it is difficult to choose \i. The algorithm is numerically 



stable, but an inappropriate choice of \i can cause 
20 instability. "In high noise conditions, the eigenvalue ratio 
of the aut6corirerat:ion i matrix : 'is : ^ ^ow* ahd' : this 'can help with 
convergence. - Ji ^ ; ; 

: t - .. . ► ; ..-*•* -v * - - 

The recursive least squares 1 1 (RLS) algorithm^ is 
computationally much more complex, but has much faster 

25 convergence than^the; LMS algorithm. It '^has two parameters, 
the forgetting factor X and the initial^'diagonal matrix teffm 
■ > 1 6. -The ^forgetting factory is set appropriate to the rate; of 
change :*of the autocorrelation of the. • ; input signal ; ' The 
diagonal" term has "little effect oil v. the- algorithm ' dhce 

30 converged, but does^affect the size -bf internal variables 
within the algorithm- during initial; -cohvergence . The *RLS 
algorithm is usually ^considered converged within a number' of 
iterations equal to: twice the filter lfehgth, which-' is 
?9:^ r 32:ly_ much faster t^_the_^ _ j^jpiitlun. ?Ji§- 
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-algorithm, can -. become numerically^ \ unstable when- the 
autocorrelation matrix of - the input .signal* is -close to being 
.a. singular matrix. • .= v 



There - are. a few - much less., common, adaptive -filter 
5 j algorithms and- the use of , these algorithms ,> has- been found 
^desirable. The. f ast.a-posteripri errpr sequential technique 
(FAEST)*' 6 algorithm . and its stabilised. version the.SFAEST 7 , 
which are also .described .in Appendix I, . have, a^.c.ojivergence 
rate close to the RLS algorithm but complexity close to the 
10 LMS algorithm. They do however, -.impose ■ <■ an, additional 
-constraint:, the input signal imus£ have<. a: siiitt* invariant 
• ; ; Property-. The shift -.invariant property, simply means- that 
■ j: t ^ ie - ipPH? ^ASP 1 *?- .tmxst- be^-tilie same; as £.he input ..signal . on the 

- - zW e Y ±oxls : ^5,?^ fpn - shif ted -on by,- - one- .sample^, with only one 
15 new sample. This property, : , is ■ not: o ;s^;t-lsJEied r by . the 

conventional architecture for a minimum mean square error 
vi. - - ( ^S E ) receiver^ fpr : _-a s ^S~fDMA u ?y^t^ A -c - ^^The, numerical 
stai?ilitir. : of rthe,, .F^ST^lsorithms inojs^^.well understood 
as for £MS .amd :;! RXS i2 b^t , i^i practice th^.;:S.EAES3? .sAgorJ-thm 
2 0 seems- , t,o . v rema.in £S£atoLe : J or , a : , suf f i cisj&t ly: riong h period of 
- 1 ime _ for. the- .^pu^pose, ^propo^ed h§r A et . s \ -The : ^ast 3»ewt on 
algorithm .(see^-. Appendix . I;) : vr-is u - ax* ; s aigpjrirthnv: - which..: can 
.sipP^fy the^ c^lculatiion^Qf : .aiiy ,tp£ ..the ^teove:- adaptive Tfilter 
algorithms if the input signal can^bjei modelled <=:£s an 
25 autoregressive filter with order less than is assumed by the 
above, filters-* - s v^-^-:^-:: =■.■ • £ 7; 

' - ■ 7 7.: ■ . , P^^ntiona^^aa^^it^ictw^: ; f£>x. the* -uplink - and 

- ; ^P^^ih^. 91...^ ^?§:T^GDMA -system-.- with, an ./adaptive ^-filter 
r. T .f 9 eL FP r : A% «b?wn^ig- ^Figjajre J,-.- r5a.-- this:, architecture,, the 

30 training r,o.f.,an^ adaptive ^FIR-f ilt^n* a ;of :leng.ehr7N-hE-l - chips 
(N being the number of chips per data bit and P the total 
. - xJ 1 ^ 11 *^ rP£ ^hipsb w 1 : qode^' ;is : donet afc:the .bit rate, using 
- ■ ^.§4§P^io?.o f#P r ; 9? - ^ f ojand . by: tho. -aigpri thm: =. to -Her: the 
. .: AJ-P* 9^nc§.^be twjeep data- ; : £E«a a« particular jtfsercf and < aaaampled 
3% estimate ^£ >6 £he data > from t£s r sou-CptXt of, the . :f i'lftez-cl; . < i . e . 
the filths feas j-a^^eff active trai^iiig:.pathtETPo.::-Th^ contents 
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■ of the filter 1 change- completely "from one iteration thereof 
to the next - This, -mean's "that ^convergence is bl'dw arid -it is 
not possible to use the FAEST or SFAEST algorithm 1 because 
the shift invariance property is not satisfied. With the 

5 IMS algorithm, Convergence is'- too slow ^and :: the time taken to 
re converge when ""a -user switches- on5 or'" oti c ±~s far too - slow. 

- . This- architecture * does work? 1 reasonably* well with the RLS 
algorithm, -although . Convergence ~ is :; Sti*i£-- hbt'-^eafy rapid and 

" tji'e computational - complexity ' is -Very -high . 3 - ; ' r " • ' : - 1 ; 

10 S\Tnmma>3ry <5T the Invention' 1 '-" ' c ; V— 1 . T...-.C* ■ .„ 

' - It- is" one r ~aitti of *thje- £ - present:^ invent ioii : to provide a 
DS -CDMA recfedve^^uSain^ an-$ adapt-ivV- : filter in'-- which"*' the 
convergence -^rs" rapid. > ir is another" aim of-* the invention to;- 
al-low the* use^ ^>f^ the - l€sk comfhdn adaptive^ algorithms'" which?* 
■15 has- not^^hith^r-to b&en pbssible" T1 : .' "-- r - r <<- ~ 

i< ' " The -present invention- provides y ar direct- Sequence code 
"division multiples-access "'"(bs-tlDMA-) ~" receiver c comprising an* 
ada'ptive ^^fil^terr^' controlled %y -an a^aptiv^ algorithm for. 
filtering data-* ^h-idh" has been multiplied by^a : spreading 1 coder 
20/ and filtered ^>y a • chaimel^ f ilter , t fie adaptive f ilter having' 
• - a length, ^apprdpriat^- to - model : %'lfe -inverse* of " the "'channel 
- * /"filter r - ; and a Multiuser 'detector :: 6peratiiig on the- output of 
the adaptive filtefx^ '^-'Z — -'jl-'^z.: ■> ^ — L 

The algorithm is preferably either trained- using' the 
25 spread-multiplied signal of a desired user only, or from a 
composite'? signals which *is^ ^he^sum^of^ tfie v ^spread-multiplied 
. r -.signals - of r^tnore thah ' on'g? - v £br ^exampJe %11V transmitting 
usersv - -This:: mearife^ that the adaptive : f il'tef wilt be : trained 
...■by - r new ■■ information at"-'*the : ' ihip -rate : 6£ r the" code! ; ' ' i: v • * 

.30 ^ In a opart! cul a*2f' etnbodiTtfeiit" of : tlie ihveht ion the ; * f ixed 
' "' r multiuser detector is :-of the"- minimum" icieari 7 squared -error 

v . (MMSE) " type r "but - £c 'rosy: -fcltie^at xvely be- of -the zercr forcing 
. ' (decorreldting) ,? c ^Volterra,- -* r Radial - —Basis • - function-, 

~ cancellat ion /".near , optimutn br~ other 1iecGdii*vj types ~ ; - 
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^ The algorithm can for .example f comprise- the least mean 
squares (LMS) or recursive -least ..squares (RLS) algorithm. 
However, because the adaptive filter of the invention 
; satisfies the- shift inyariance. -property, the algorithm may 
5 alternatively comprise the f ast .^posteriori error sequential 
•technique (FAEST) ..algorithm, : . -the stabilised FAEST (SFAEST) 
algorithm, and the above algorithms or .others may be used in 
combination with the Fast Newton algorithm. 

- Brief Description of the Drawings rcv *= ~ 

10 ; , _ In v order that- the present ; invention .may be -more 
readily understood,, -reference will now .be^made,. by way of 
example only, to the accompanying drawings, in which: - 

Figure 1 shows the conventional- DS-GDMA . architecture 
already discussed ; 

15-. - Figure .2. shows DS-CDMA - architecture according. to„, an 

- : embodiment of t'he_ invention ^ : L ^ * - ^ >■ - »• ■ 

Figures 3, 4 and 5 are graphs of simulation results, 
showing -respectiviely r; the : comparative, -convergence of the 
architectures, - of , Figures 1 and 2, the relative, .convergence 
2 0 rates of different algorithms using the architecture of 
Figure 2> and :7 the-- bit .error : : ra : tio : (BER) ; results for the 
architecture of Figure : 2-;.— : v ^ ~ z _ c = ^ _ 

.,3 Figure : .6^ schema : t^ system with no 

channel model; 

2;5;- -^:;r Figures^ -and.^^ are ^graphs ,. showing signal- to noise 
• - ; ; .performance - of:.. i aa^iener filter Ga r lculated for 7 chip : and 31 
: ; chip .-.Gold-; codes _ respectively;; w : - : iS r m , • ^ v 

; ■ .: Figures; r? ; and 10. are graphs -^showing .- bit error r rate 
i3 ,(BE;R) per^qmance of .the r Wiener .filter calculated. fo?r .7^ chip 
3^0* and- -31;. chi© Gold -codes -respectively^ - f - T . r :j : - ;t : 
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Tigure 11 is a graph* Showing convergence properties of 
the LMS and RLS adaptive "ff Iters ; - * - ■' 

Figure 12 is a graph" showing BER- performance- r 6-f- the 
LMS 'and* RMS : ' algorithms ^after- 1000 * iterations allowed for 
5 Convergence/" compared- "with the' Wiener* optimal and -a matched 
filter with no- MAP ; ° - * - v ~ :: * :c r ' '' ■ - " 

Figures 13 and 14 are graphs plotting BER against 
number of active users in- ah -AWGN 1 channel* ^alrid a - -stationary 
■ multipath* channel - respectively, iall users being equal power 
10 and the" spreading code'' length being : 7; " ' - .-.--v - .1- 

Figure 15 schematically shows the construction of a 
' received' sighal' Y(n) ; " ; - J ^ ;: ^ 7 ° fr:.vr "-v ,y 

Figures 16 a) to d) schematically show the structures 
of a"- matched filt^, ; -a parail el 1 canceller -using matched-* 
15 filters, a Wiener filter arid' "a'^psLrall^ei"' canceller -'using 
Wiener filters respectively; 

r: Figiire ' i-7 ::; ' is- a - graph showing -the BER <peirf oritiahce of* 
- :: the - : f ilter^'-silown iri-Tigures* 16 a}"" ; to d) ; and * — - ~ : - 

Figures 18 -a:)' £o : dV : aire * : gf aphs plotting simulated^ BER 
20 averaged over all users against- number of users ' for-" -four 
differnt signal to additive Gaussian noise ratios, with 
- : 6-6 ; : 00b -data bits' per • user arid* a'-seqSen-ce- length of 64. 

Deatiled Description of the Prefer ^^ ^odiments 

25~-r-r. cz pi^r^ 2 * : shows -^-DS-CD^ - ttari'smitter- -&nd - receiver 
- c architecture' • corr^rd^&ing" spreading- 1 * 1 mfeahs" -'2^ respectively 
operated by each user in which 'a "ddta"'* sighal :> from- the": user 
is multiplied by one of a set of spreading codes uniquely 
~ ~ allocated "to "^th^ ; 'Respective userP : $fie dat^-i-s. supplied to 
30 "the spreaLdihg'means-afc- its bit irate 'and fefie ^ddde : 'is input at 
its chip rate, therfc^beirig N : 'code chips ^£&£ ea'dh da^ta-bit." 
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. ..' , . The spread-multiplied data is ^ summed at 4 and filtered 
. in - a .fini.te .impulse response .. (£IR) .channel filter 6 of 
.. length P ..chips, the output .of which. „is transmitted, with 
incidental- Gaussian noise/. to. a .\ receiver comprising an 
5 adaptive . FIR .filter . 8, also, of., length P chips,. 

Ah adaptive, algorithm ,10, controls .and trains., the 
_ adaptive filter 8. The algorithm 10 basically fi n ds an 

- error signal, which is. the. difference between (a), either the 

- spread-multiplied data from the. desired , user _only or the 
10 composite ^spre&jd-mjiltiplied signal., and (b) „ the . output of 

th^., adaptive r f ilter 6 . The., choice... of training . data in (a) 
can either be preset cr switchable.. .The. .effective training 
paths ETP for the adaptive filter, representing this 
„ operation of the algorithm are shown in figure 2 . . 

15 • r . ... . The ., adapt iy©.. .filter . . 8 r _acts. similarly to~ ; (but . not 

..exactly the,,. same v „ : asL ; a cpnventional equaliser^. , As 
. previously- ^ stated training. ; can : b,e .pn _the .desired^user's 
: signal only ,. or„ oji the .composite. ..chip., .rate, signa.1 9 , The 
. f ilter .8 ..has. ; less/..tap§.^. than the ..adaptive, filter . l in., the 
20 convent ional. v .archit^cti^r,e .and,, is trained -at . the... chip rate r 
which, means much faster convergence ..The- adaptive ^filter 8 
als.a satisfies,, the shi^t,, inyarianc.% property ^required .for 
the FAEST and SFAEST algorithm and this allows „ _LMS 
complexity levels with RLS convergence rates . A fixed 
25 mul,tijiser. d^tector^ i 14 ...is precalculated and stpred at the 
receiver f and. is r: bas^4 :£ pa,topwledge .^obtained about, the number 
: : of ;J .users £ in ttie . system and-, .their r spreading codes.. . _ This 
h : detec£pr. l ; 4-.can r b.e o£ xixe. MMS-E type 9 (used., in* the comparisons 
......here)- or it,. r ^carube. .<pf.^ a. different,. .type,,; for example zero 

3.0. .forcing .or .decorrelating 1 ?; -Voterra 11 ' JRadial Basis function 11 . 
r . cancellation J)ase4"';. " r pr, n^ean. optimum decoding, .based on 
. ... -Viterbi- decoding 15 . _ ... •. - ~ . _ ... ; . ... . > 



- j- Suanmary o£- Multiuser Detection Techniques 

^ T \ ts - r Til ^i ^convent i onal js ingle, user, det ec tor or, spreading 
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code matched f ilt'er " cafcuiates" the correlation' between the 
received signal kiid tllef spreading code (or the spreading 
code convolved with the channel impulse -response *ih~a 
imiltipath channel) 6ve2? t*h6 , data' bit period. The 'multiple 
5 access interference ''CMSlf impact's: the' ability to recover the 
desired transmissions with this receiver and is more 
pronounced for high power interfering users./* " 

"The matefied'f inter : "i r s'" t:hk : optimtim detector in AWGN : but 
the' MAI ' crdb*b-correi*atioil* : 'term^' acre' not* Gaussian*'' unless 
10 there' kre 'many active '-us 6*3" r Irr t r hiS A sit&ktitrfi the optimal 
detector is" Hot' the" conventional matched "filter** but rather 
" a form of r 'ftiuiti-user deteetdr;- 1 '' ' " " :r * " * ' * ' :: 

The* d^^btor : th^t? r yielc$s' the' "md3t likely transmitted 
sequence maximizes the probability that it was transmitted. 
15 Whefi all ^OBSib^e" , transmittre-d 7 ***seqTifehces'* " are 1 equally- 
prob'ab r le, thil is^fihe maximum f ike^hodd r ^egu^nce e'stimator 
" (MUSE) . ~fth£n "^ftrpf emeiited : with " Vitefbi ** ^algorithm J the 
complexity* iV%5cponentiai- i^- rhe nfimber of° users r The MLSE 
■* 'must also 'VsriTt^te" th^ it 
20 ^ foWrs^^tlle r mo6il'e' "tis er transtnittf^r : p6wet J ' fcontroi accuracy 
r fe<^ir^ment Despite t Se" peirfbrmhfe^ and capacity gains 
* ; : o^e^" : ddnv^ <iet : ecribn:, '"the 1 "MLSE is* riot ; a 'practical 

; " •'"soiutiarf. " "' Lz cr "* : ' ' : " Tl * * lt 

:rs r He£i" th€ irlsxcr -forcing ^dr " ^ecdrtelating i J detector- 
i2S^irtpl^ment!s 0 £ri inverse*' of -tfie v ^cei^cf ^ii&nar aufocorrelation 
' ■ matri* ei (n^ -'Mrdh-'- tfiat the J output - is 

* ' completely^ J f ree : * of multripl^^ adhWs 'interference s : *The 
decorf§l ; atirig- detei'ctor~V^ in the *"i l ate 

'1970s* and it provides substantial '^perf orm^ce?capacity r gains 
3 0 &ver°ttie : conventional ^&etrector - and doesr not need z to- knbw in 
advance the received signal amplitudes/ avoiding" sensitivity 
to estimation error. It has computational complexity 
s igni f i c ant ly lKWe ±*' -MrMLSE e >r ^ofch^ i s 

that it corresponds to the MLSE when the user energies are 
~3 $• VrSdSSwnT 0 A disa^rkntage of 3fhis detractor v rs°th^rt r it causes 



WO 00/51260 



FCT/GBOO/00649 



- 8 - 

: noise enhancement i.e._ the power .associated with, the noise 
..term at. the output of .the decqrrelating detector is always 
greater or equal to .tlj.e noise .term at the output of the 
. conventional., detector. A, more, significant, disadvantage of 
5 the decbrrelating. detector is .that, .the .computations needed 
. to invert Xhe matrix are, considerable . 

_ -. _^A\pp.ssibly_..^ ^eceiyer^stnicture is to 

build an adaptive filter, that; minimises (at;, its. output) the 
error, power* \This implements a partial, or modified^ inverse . 
10 of the autocorrelation matrix, dependent on the, level of 
background noise, to balance the desire to decouple the 
r . users for MAI reduction _while. N no r t enhancing the ^ background 
noise. Again this receiver structure . implements a matrix 
inversion operation and here we can apply the recursive 
15 .-adaptiye_f.ilter .technique^ T^is^ detect ox. differs from the 
..decorrelating, detpQtQr^.in that.^i.t takes the^noise terms, into 
account to minimise the-^pise enhancement .... „Thus if noise . is 
low the minimum mean squared error (MMSE) receiver 
approaches the decorrelating detector and achieves, qlose to 
20 optimum performance. On the other hand, if the multiple 
.access interference is.smaH .compared to, the .noise,, then the 
.. matched filter.. detectpr, solution .is, apprpa.che.d- This is 
exactly, .analogous to .the . MMSE. linear, -equ.alizer used, to 
combat inter- symbol. int,e3rf ^^.ence.. Uxilike..the de.cor^lating 
25 detector f . it require. s. ^sXimation. of , the . received, amplitudes 
but -its complexity, is._indep.end^nt .xjf jnumber of active users 
. .and .explicit, knowledge., .of -the, CDMA ~ spreading .sequenqes. is 
. H .not required., for an ^daptive^ implementation t _ . - . 

3 0 Other researchers. have. rt investigated., radial, bf. sis 

function (which in its complete implementation is similar to 
. .. MLSE,-.. various. -s ; implif ications^ have .been., suggested) and 
. t Vol.terra. nonlinear, .approaches -.(where ,.the received sigpai is 
. .passed. through. .-a. jp.ow.er. .series . nonlinear expansion, .before 
3 5 . applying : a linear, filter :..eg*- decorrelating. or . MMSE) . ^. Note 
. -.that if a .MMSE, -error. ..type .fixed multiuser detector is. used, 
s the-, .overall impulse response of -.the ...adaptive „f il.ter followed 
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by the fixed detector is" hot 1 ekactly the ' samfe~" a's"'' f or the 
conventional "architecture "even When both systems are fully 
convergecL This will gfive different bit ferror * ratio "(BER) 
results " even whetrT both 1 'the " adaptive filters" in the" 'two 
5 receivers are f uliy 1 "converged. 'Note iilsb th^t the adaptive 
filter 8 has an extremely "Tow signal to 'noise *rat*xb at its 
output as the processing gain has not been applied at this 
stage. ** This' leads to "a Ibw' eigenvalue" £pr£Sd which can help 
the -convergence' aiid~ Stability il bf Tbdme a'dkptive algorithms 
10 but others' ei&ibit 'iri'st'ability iii the extremely high 'noise 
conditions' ""^ ~ " w ' :; * '* . 

Posisl^le^'multltis^r v ^det:'ecrt:orj^ ' are~ 'discussed'* fn* lJ mbre 
/detail iii ^ppeftdix'Ti; ' "'""^ 1 L ' — — - ^-u*^ j - "c ; - - r- , 

~ Thfe" blitputr ; *£r%Trf the dstefctb£ 0 i4 :t is" ; *d^ 
15' bit 1 'rate ^ at : 't&W synchronous" pbiht?^~ J tS^ive "the '" required' : 
estirtiatTe of r the" uWer >: s ' datm M'grikrv' " bc ~ • - ~ — r - ~ :Zi *- 

S imulation Results ** ' ' * J ~ ' ' * •'-* 

Ail "slmiii at ibn~ 're suit's p3rfes'In^i&cf r bfeiDW a are for^ie chip 
spreading 1 idotife-S"'" (i :e . ~ ' P=i6 j * : arid' ~ %: J six tap stationary 
20 ch^ifiher:^ Firstly* Wer show %hkt "fchie hew afbiiitecture has far 
"f ££ter°bf c driver ^eiice th%n~*th;£ ; lidii^ention^l" one : . Figure 3 
' ' shows ^ cotivergence ^curvfes^ f or Vher ccftivientiOTial architecture 
: * 6f Figur^ i^an&^the ^w^axb^tebHix^ 
. architecture 1 i^^ftiuc^'^f^strel: "tkaii-^ thl or idbnv r entibrikl : "one, 
25 mostly becaus^'it'" ^ of" the 

bit rate but also because the eigenvalue ratio of the 
aAtobbr relation* 'iV v "re v ducecf "in"* hi^ir J rib&ser ? : rr ~- - w 

\ : ~~ ^Th& graph ""of ^igtire^ 3' shbws 7 eri&^ttole c? Averaged squared 

eari%r tlie' "bu^^^ r ~/bf 

.30" architectures plotted 0 ^ chips; *for 

the single "user * -censev 1 CD So th* r 'curves" us^-the *LMS "algorithm 
* * ' VitSf tfhe ; 'value of "7? ^ttdfvld^lly- :; o^€xmrsed ^of each" " : 7 The 
" r '" ii'sW "of 3: chip tate^trairiing ^nst^ad^ bf - : bi't r rate"- training, 
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. increas.es the convergence, by a factor of .16, but in fact the 
. convergence is more than 16. times faster because high noise 
reduces the eigenvalue ratio in the autocorrelation. The 
. new algorithm converges to a MMSE 1.6 times higher., than the 
5. conventional, architecture, but this, loss .will be largely 
regained through . the f i}?e4 -multiuser . detector . 

_ ..Figure 4 shows, the relative convergence rates of so me 

different adaptive filter algorithms using the architecture 
of Figure 2 . The use of the SFAEST algorithm is only 
10 possible in the new architecture, as the conventional one 
does not satisfy the shift invariance property. 

Relatively ^speaking, the LMS algorithm is very slow. 
It can be made faster by increasing the value of \i but when 
this value is increased too far, the LMS algorithm does not 
15 converge to the MMSE error floor (mis adjustment error) . The 
RLS is fast, but at the price of high computational com- 
plexity. The SFAEST algorithm is as quick as the RLS 
algorithm if it is initialised appropriately. 

Figure 5 shows BER results for the new architecture 
20 after allowing 160 chips (only 10 databits) for convergence 
of the adaptive filter. 

With only 160 databits for convergence, the LMS 
algorithm is not fully converged when the training period 
ends and this results in a slightly poorer performance than 
25 is obtainable with the conventional architecture. The RLS 
and SFAEST give similar results. 

The architecture of the receiver of the invention is 
much faster converging and tracking than the conventional 
architecture with only a small loss in BER performance for 
3 0 a stationary channel. In a time varying channel, average 
BER will be better for the new architecture because of 
reduced tracking errors . The fast adaptive algorithms 
(FAEST, SFAEST, preferably in combination with Fast Newton) 
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allow the new architecture to achieve ^good convergence 
without' the computational' Complexity of the RLS algorithm. 

The receiver or the' invention cfould be integrated in 
a "hard-wired" form or could be made capable of being 
5 updated by usincj' fecorif igu£abl£ ^br replaceable firmware:" 
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APPENDIX I 

; r Summaries of - .Adaptive' Alcroa-i t-Vim^ . 

• Least mean square algorithm '• . '.: ■ • . ? 

- ■ . From, the stochastic gradient approach-' [1] one obtains 

5. the well known least: mean squares (LMS) --algorithm which was 
first s introduced in. ..i960 ' by Widrow arid Hof f [2]. This 
algorithm is uncomplicated' and yields acceptable performance 
-= -. in most -.cases. .. The LMS - . technique "' is : probably'- the~mosT^ 
• ''frequently used adaptive algorithm in current communication 
10 syst.e-m.. It is derived -by applying -the method- of' steepest 
«' -i descent •..minimisation •:. to", the. Wiener-Hopf - equations which 
define- the -optimum- -Wiener filter- and one' -"obtains a 'simple 
. -.recursive, scheme of .it he - f orm -[3 ] ■ ' ' '- ";-:.'-' -.. '•; • 

: ; J i ;.; ; ne * • _ W / - i?ld - • ^ - \ + I inning Inew .] 
'l ! ■ 1 W W^- /., - ■} ^ ^S^ /■ ', I . Iff J I iAformatipn J (1) 

where the new information, consists • of'^ t-he -"product of' the 
15 filter input vector and the error signal, i.e. the 
difference between the desired: falt'er~6ut^ 
. output. :of : the :fiater-.s^sen-efally'-6rie-- 'can? express the LMS 

....algorithm (and 'many mother adaptive ; -alg6ri-'thms') by means of 
■;. ■; three, expressions- C'l]-:~ " - • ~r.-6 z ; 2..;- 

20 Filter- output : w(t) T x > f) ; ;. ■ ,: ': ^ 

Adaption ""error = • ; ^) 

' - - -> - ' ■ '-"."^ Tap:, weight : update :-w(fT^= vfcK-r /ix(l)e{t) ' ~~ - r ~- ( 4 ) 
where x ;is rthe:' ifil^e^r^nput' vecterf- w r the- : -tx^n-sposed- tap 
, ..; weight ..vector, .. d ^represents- the- 1 -de-sired-- filter- output -'and fi 
25.- is the; iea'rning-: rate ('sl^e = p"^i2e-'param¥t : err.-- ; v.^- " 



■■ ■" .Several- v modificati'oiis 1 : "have'' been "made- to improve the 
: : -IiMS' xalgoritlhrft', -'- -the- "most;- ^ important-- - one' • resulted'- ' in ; " the 
a. nbrmalisedvEMS :*-NLMS: which normal ifefes '-the -"'adapt ion 

.: error e :. u'singv instant ahebus^ estimates'- of 1 the input '-vector 
•30. power - ||x| •?'..-.: Doing e so • -largely -improves- the ' algorithm' s 
r r. . stability :cfc3racce;ris£ics -and -ailo'ws faster -"convergence :" The 
LMS algorithm is computationally simple but its convergence 
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rate is slow ^nd higl^l^^depe^dent. ;on the. properties of the 
input signal, more specifically on the eigenvalue ratio of 
the autocorrelation matrix. When rmany elements -of : the* input 
signal are* unknown-,^ for .'example the - channel .in' a mobile 
5 communications., system, - --it is diff icuit." to; choose- The 
algorithm is numerically stable? -but. an inappropriate choice 
• of ^ can, cause, -instability, ^Tn- high' :noxse conditions., * the 
eigenvalue ratio ; of ,:.the;;:autocorr^lation matrix /is. low and 
this --can help-.- with :: _ -convergence . .Ottier; r versions j of the-- LMS 
10 such as leaky K .~ K signed;. ,-anci quantised LMS1. [3f, were, developed 
but % only the - -NLMS.v. achieved ^ * sufficient ./stability at 
acceptable : ; perf orman-ce-. : arid ' is. itheref ore:: eon-side'red ^as a 
candidate for application - irir advanced r.muitriu.ser -detection . 

"To illustrate*"-^ the_ ^improved . perf ormahce . of the 
15 normalised "LMS ^ as 'compared tor "the stahdaird LMS' we will 
present the performance of both the LMS and NLMS. later in 
this • document^ ^report r. • .10 L •j^.irr : W; : ;.. y. j y:^r.- 

. Recursive least square algoasit-iro; =, :: 1 r.^evc = 7 ' ■'■ . ■ 
. % . ■ , .-The^ : mpst^ r pq ^technique; is is 

20 prpbably the^- - : r : ecursive ,le,ast square s\U RLS algorithm, The 
RLS algorithm is computationally much more -complex, but has 
much faster convergence, tha^,. ther.J»MS.c algorithm. It has two 
parameters, th.e forgetting- factor h&yanjfl the initialisation 
:-i factor -6; t pr-thej diagonal ^a-tr^ see" Equations below. 

%5 : Th.e, _ ^rge ; tvting f factor set- -4PEi^pjriate:-tOi.:-the ^ rate of 
change ; v of _ -,tbe h autqc^rrelat ioni^Qf ~^ct*Le o inputs -signal ; : The 
diagonal term : .has ^little, ^ef £gc£) &ns rXhe algorithm i once 
converged, but does effect the size of internal variables 
within r; : the ; aljgQ-ri.ti>jn : _ during-; ; ~kn^^^^ The RLS 

30.. algorithi^.is usual ly : -considered: cpnverged .wi^tttin ^number of 
. , 4^;er.ations *.egual rfr to twice the*' filter clettgfch*. .:. JwrhichT is 
. : generally much . f aste-r^ -than vet he- r hH$^.ral^£x£hm. .The. =■ RLS 
. ^.^gorithm carV; ;? become ./numeoricaliy ;r:unstabl'p; 1 when- the 
, j 7 aiitaqo^ :inpuft v.Signal/:i^i:close' ±o* being 

35. _a..aing}ilar ^matrix . ■ r - ; r ^. •r;^;;;':Xw z j: rr t ?i ?xr.o~~ '.z :--\ f - r 
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We can summarise the„ standard RLS .algorithm as follows 

.._.„ ._ ..... .. .... _ . INITIALISATION - - 

w (o) = o .;. 

"• — SUBSEQUENT" ITERATIONS' ' 

: . . ... . Z^z^r^t^)* r '■ . . ■ n 

' • - w - '• P-00 -X (?-(*-= 1) - S (<) !))• '■" -., - : - ' (8) 

5 _ _ T he - ve . ctor . ? i-- s - - I S now i?l a §-. .the gain vector; of the 
. ..algorithm, .Pit the - inverse-of ■'- thV correlation 
_ matrix _ _ _ „L 

which can be computed in a recursive .manner' 6y ; 
„ ; -M ai i?,, 9? Jr£?.*^^ Aemrna v^Th^.j^r?iable. e,: denotes 

10_tJi:e .a zprdLori '.-error, of -tttdi'i il ter-es t iniaii^on-and 'the 1 filter 
weights are again represented by the vector w. 

'" "Fast a -posterior i e^or sequential* technique (FAESTO - \ 

..^e ^A®*\ ,^.rP°?^^ rio f error sequential* • technique 
(FAEST) has * first, been reported, by Carayannis- et ;4 : al- : [-5] ; in 
15 the.' context" of least-square (LS) filtering. .Compared to 
standard LS algorithms the FAEST uses a different approach 
for calculating the Kalman gain vector. p^sed\ Cx &vr: : the 
a-posteriori_ er?gr. ?qrmulation . ^ rather ^ -than Jthe a-priori 
" ^error ^formulation as .iJJsejd g £p f ast -.Kalman algorithms./ [6] . 
?A Assuming the^ shift inyaf i££ce ^property of z the input signal 
J- and ^ in t r pducing ^ a slightly 2 m6di£ ied r -vers ion, of -the*, Kalman 
S^in, ^he ?AEST algorithm.^ -manages, to. perform a.-: direct 
updating of the kalman gain without invoking matrix-vector 
multiplications. 

25 The Kalman gain update according to the FAEST 

algorithm is summarised in the Table 1. 
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Table i : ' The J FAEST algorithm." 



Time update of gain vector 






Multiplications 


Divisions 




(t) - a: v (f - 1) x,v (r- 






•V 
















1 


a' v (I) = Xair (t - 1) + e' v 


(t)si. 


(*) 




2 






rA'rl (*) = 




/.v out)- 








2 














1 




WA>1 (f) = 




0 


e 


' v (*j • 


■ l 












w.v (t - 1) 








1 


a.v (t) = a.v 


(« - 1) - c'v («) WAT (t - 1) 




1 .v 


- 


e'v (*) = ~ 




I? (?) - 






2 


- 














1 




r ; v(t)r«(t)mi.(«) . 










L- 










- 


<-. :.. ih -- 


a s v (:) = Aft» v (t - 1) -r e? v 


I*) *°v 


(*) 




2 






0 






-b N {t-l) 1 


iV .. 






a.v (i) = b.v 


(i r - : l) 7 4(:)w^.. 


(«) 


..... •■ :.-=■■' "J: ■' 


Z.VH.r iV;.- - 




Total numbe 


r of multiplications and divisions , -a.., 




,5 . 



It is well known that the FAEST algorithm^ suffers, from 
seve-ir- stability -pVbfciyws and we wilF' there fore' not describe- 
•:• - this" algoxi'thm : -in more -detail "but 'focus" "on its ' stabilised 
5 Version-;-- the^SFAEST . For mdre" info rmafc ion about FAEST please 
r-ef -er -'to [S-]-!" " '•'-■■ •• r-.-ps-Jii*. 



-' Stabilised FAEST ' - r..... . . _• . .. 

1 ^ e &t : aHilised- r version- of the FA^sf "algorithm" has been 

' : c : derived - fcy- [¥F '-and ' is " : preseWed w '}rable *"' 3 ". " For "easier 
10. reading- of" the ' equat ions buiidln§ ~€&£ * s'tabills^ed" ' FAEST 
• ? - • S FAEST )- ' ; a-rgbr'i'thm; - ' we^°^ii^ r ;fest " listf *' ail" "variables 

• . occurring arid explain; ^n¥ir~meafilfig fi wiVh"/ a .few words','" see 

• "Table --2: ; - : <-"-'- •- * -' : " 
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Table 2: Variable def initions :of . the -SFAEST-- algorithm. 



Variable Name 


- - ^ Definition - • 




... .:- - Filter input-at time-t- • 


y{t) 


Filter output, y(t) = OVf < 0 * ' 


x,v(t) 


Input vector [x(t), . . . , x(t -» JV" + l)j : 


- 


* - *■ * Dual Kalman gain vector - 




Forward predictor coefficients 


— - 'biv(t) 


rBackward predictor coefficients- y 


h,v(t) 


" ""Filter, coefficients 


• e.v(t) 


A- priori .filter err<?r ; , , 




A-posteriori -filter error 


' ev(t),eirlt) 


Forwar a predict bx "error anH its corrected Version"' 




Backward predictor error and -its 'corrected version 


aV(t) 


Forward prediction error power 


* a*v(t) 


; - / ^ -Backward prediction: error, power ~ . i ' . 


7iV(t)^7iV+l (*) 


0 < 7 ,v(t) < 1. . ■■ - 






"tvW^.vtn - 


. Difference of. errors. atuT the correctecCerrox difference.. 


A 


exp. forgetting factor 0 < A <T - 







Table 3 shows the successive steps : " 6f\ J the SFAEST 
algorithm and also evaluates, the * . cqmplexity . in\ teirms . of 
multiplications and divisions required:-; . \: i \ 
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Table;*3':A Tasks of . SFAEST . in order of -computations 
and number of required MUL/DIV operations 



*.=.•.- -.- .-- -Tasks -of SFAEST- * , • , 


MUL- 


DIV 


Available at time fc- ,~"T.. — : „• - • - 

&tf{t - 1), a*(t~ Di'KirXt—.i), Mi~i) r *#(t---- l)r - - 

7;V (t - 1), a&<*- iM&(*- 1)-— - -' ' - "~ 






New Information:---- « — . -~ L*~. .... . ____ 

x(t). v(t) ._. : '.;..J. ; .t.±."-" • J' . ~ . ..... . : 






Computation -of the residuals and,.corrections: . _ . .. . . . _ 

A-prioEirJ*orwaxd/.iax±w2^d ^. ^ 

1) eL{t)~= x.(t) « a.v^t - r).^vX^ 1)1. : _ J I 

2) ef v (tV.= xCt -3)^Mtf _ [\... 


> 

• n' 




Normalisation T parameters :;- : I : f ^ '-"'J 

4) 9.v"L<) = i + 7/v^ft)4Msg(^j) -i^iy^vji : - ij?J _ 

5) *wft) -a-""?-?*** aa4-?v(t>-^W/!\'-- / > ■ -'•-»• 

6) Jb:v,(r-:T) =-A-! v 7.v(^ri)a^(t =-1)- _ : -•• -.--•=r-, • — ; • 


'.. 2.- 


1 
1 

5 


Difference "of gftots?. ;,, -7-;-- 

7) £v(t) ^^-^(f'-^'W^^'ft -• : i>g^(t= l>~ - ~ 


— 3-- 




Difference of errors of correcred filters 

8) £v(*) - i-j-pd— r jV (t))-r^V(t-X)(l-Tr.v(t-l)) 


4 


1 


Corrected a-priori forward/>ackward .errors: T _ .... 

9) • e£(*j - (1 U w(* - l))fcv(t - l)p^(i) "~ ' " ; 

10) - ^(t) -= : Aaf v (V^ f)-7-7^0 - 1) (^r(t)')" 2 ^ * J " * " " ' ; - 

11) 4f(t) m ef v (t) - (1 - 7i*(*)*P?w(<) ■ ^'-" v - ~ — -- : 

12) a» v (t) = Aa* v (« - 1) 4- tv(«)#v(*)) 2 


' ' "3 . 

' 'i : • ' 

- ■ "2' r : . 
3 


- 


Time update of the duai Kalman gain g;v(t): 
Extended Kalman gain vector: 

13) s,v + iW - | g^.D " Ao / v(t . :) -a, v (t-l) 


N + l 


i 


Forward filter: 

14) a,v(t) = a/v(t - 1) - 7,v(« - + M* - l)/£v (*))&v(t - 1) 


iV-i-4 




Dual Kalman gain: 

is) [^^s^w-s^wf-^r 15 " 


iV 




Backward filter: 

16) h s {t) = b. v (t - 1) - A*(«)(e 6 w (t) + P? ( v(0)g 1 v(t) 


iV-s-3 




Time update of the filter hjv(t): 
13) e jV («) = y{t) - h%{t)x N {t) 

19) £Ar(t) = 7<v(*Kv(t) 

20) hv(t) = h,v(* - 1) - «/v(*)giv(t) 


1 

N 




Total number of multiplications and divisions: | 8JV -i- 36 
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Initialisation, stabilisation and optimisation issues 
related to the SFAEST 

The value of the forgetting factor X largely depends 
on the rate by which the input ' signal changes and defines 
5 _the„ meiiKjry -length, of the- algorithm.. The eigenvalue spread of 

. the weighted ..input covalence matrix given by. .- - - 

' - - ; " "t 

\ " . tzzQ ^ ■ 

plays an important part in.. determining k ■ A . suitable value 
could for example be A = 0.98 but the type of the input 
signal ...needs tp be considered; ^and : it-;is most likely that 
10 . changing, J^e^yreen stat ic and- dynamic users-,, and v static and 
dynamic channels will require -different values -for . X 

_ c . . . The .chQice,. fjcx.the, variable p ; Is, not -straight forward 

P-°^ . r }^ e -}, . f 9^., 9' alculatin ^ an appropriate \ : value 

exists^ In, L7 : : : i,t'r. is-jmentioned^ that- a initialisation; .can be 
15 done by mear.^ of, :> an .estimate such , as \- h L : r.- : -: :jb 

h .1 2 i .: with . p 0 2=.0.Q5,: ; ..^ '(10) 

This ; cpu^d however not., be., confirmed ^y;,simu; at ions- v -carried 
..put ..in th^ .course , of , thi ; ,s work.. Usually -a value, of : //^=t^ has 
. . .feegn chosen .but an adaptive; est imatipn : -Qf_/> .during pp^rat ion 
of the algorithm £ might : pr : ove ; •advantageous;... The most; ctucial 
: 20 .parameters were., _ff ound 6 tOr be- the - f orward/backwarqt ; .error 
_ powers.. Initialisatipn :I has.- -been _.:done--jay means:c ;Ofr the 
following rul.e^ ~ ^-^^..v^c" ~ ■' * -'-^r \zlz 

*^(b^^ *---^ a - s ~ ■' : (i 2 ) 

As ...qan . : be.^seen 'both ejror powers depend on the . value:.. /x and 
^ en ^^: ^ is ^ .^ a ^ m ^ e ? ed §.z t 9-. be- initialised with, ..care. A 
.25 typical., ..stable ;J range : f*of..: yas j: f oundL--to.; be ,betwe.en 10 < \i 
< 100. From a stability point of view it is ^advantageous, to 
monitor the evolution of the variable 7„(t). To prevent 
. : ... -diy^rgenoe £ af^ 7 the .-.algorithm r , r s this- ^variable. ; should be 
, . r^str ? i^t^d :j tp values between and-i?;; A ^ecottraended, nml-e [7] 



005126QA1 I > 



WO 00/51260 



PCT/GBOO/00649 



" 19 - Appendix I 

'to assure convergence' is r; £o reinitialise "the algorithm' when 
the condition ' w 

- i ^ :.i>r^ &&>v<&V) , - (13) 

becomes" true,, where v is™ a" 'small constant : "This precaution 
will also take "'bare of- the case where "the gain tends towards 
5 zero, hence 7^ -* 0. -x*"'' 7 .>:-S\ 

Fast Newtba transversal filter v ' ; " ^ 

: --"'-The Fast Newton*" al-gorifchfrr is 'kii algorithm "^hicfi can 
"■ simplify the calculation; c of : any of : the sboVe adaptive filter 
algorithms - if" the-' "input y signal can be r m6deired* as - an 
10 autoregressive filter with order less than is assumed by the 
above : filters : '- Fast- Newtorr' transversal' ~fiiteirs originate 
from the area of speech enhancement" atnd echo cancel i'kti&n. 

• • Their-' -main- -feature- is- a : fast caiculatioh of the gkin ^IcHor 

as required in many^S' ada^iVe' : Sgorit^s'. J£,r " • - c - 1 - 

15 The stabilised-fast Nevton;;t"ransversal filter (SFNTF) 

is essentially a computational accelerator for any least 

• - "square - \^LS Y- -algorithm . r *~ ' Operating : as" "k~* "higher- level " 

-adaptive -"^predictor ~ -it can * \is*e any t; LS* " algorithm as a 
~ subroutine' within its-' own^ algbrithtrf. r However, the order of 
20.:this-X.S' filter- -Sari 'be 'chosen :: 6a be^smaller * : than EhS aetual 
r filter " -ord£r- - o^ ^ the ? SFNTF- ' J afgorithmT 3 ?A : sophisticated 
• ■ predictor part ^fchen -^extr^blaiies* ' £ fche' ? "remaining' '"filter 
coefficients to gain the complete set of ^litrer 1 Coefficients 
. as required according to th| definition of the SFNTF length. 
25 it is this feature that should, : igake the SFNTF a potentially 
attractive adaptive algorithm for many application and the 
v: usefulness' * of- SFNTFs- fdr'adv^ be 
. examined in -the -f littered As :: f Sir- is tEiV-'^epcJrt"' IsvcSnttetned, 
•we will- however resT^ the 
30 -al^o^ithfn only r • '-^h/.'r:: -:i - i-l V--" - ^ £ r.^ : . 01 : 

ilucm Tor discuss the-deta21s j: 6f ^the STW£F ~v?e w±tr k^axh^'list 
aftl- va^iat>les- -as pre's^aStr- ±h : &Te : algorltlurf* wftlf Vom^irbtas of 
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explanation about the meaning .of the variables. Table 4 



lists - these definitions 



Table" 4-: ^Variable definitions of the SFNTF ' algorithm. 



' - • ■- '■■ . 




" Variable Name j 


Definition - - > v. 


. ,-.<pfeK- z.\ 


-Filter input at time.i; * . - ; - : 


y(tl , . . ! 


Filter output,. y(t) = .0 Vt < 0, 


x : v(t) " I 


Input vector [x(t) t ...,x(t - N + 1] . 


" S;V(t) •' ' 


"Du^K^hiati''gaffi^ctor of" order jV" 


qjv(*)iqAT+i(*) 


temporary vectors to compute g,v(*) in version 
2 and 31 - \ " 




temporary vectors to compute g/v(t)- 


a.v(i) 


Forward predictor coefficients 


b.v(0 


Backward predicxor coefficients 


h,v(t) 


Filter coefficients 


R-v 


Input sample covarian'ce: matrix 


r.v 


Vector that. builds R ; v, see end of tabie. 


e ; v(«) 


A-priori filter error 


•: e 3 v («) ...-' •: 


7 At posteriori filter error.;. ' : . ' f I 




Forward predictor error and its corrected^version-- 


<v(*),*?v(*) 


Backward predictor error and its corrected ver- 




sion.. , 1 


a'vW 


Forward prediction error power 


*°v(*) 


Backward prediction error powerr 


7JV (t) ,7iV+i (t) . 


0<lN(t)<l 






€iv(*),^v(<) 


Difference of errors and the corrected error dif- 
ference^ - r ;~ : ; - .. ,:;-..::v .:. vr - A 


A 


exD.^forgetting- factor 0 <„A <J_ _ . 



We can now describe the' "S?NTF' : with "all its ' equations . 
5 For easier reading of the following equations, we first 
define "two neV~^e*bto"rs c s : "knd 'u: " c *" 



' +i(t)w A«j,(t-i)l i r. 



(15) 



Comparing th6se.-def jLnitioTis v wi-th Equation 4, we notice 
the similarity with the" updating--: part ^xe. While A 
represents the f orgettihg" "factor similar to' z 



the second 



! - : 10 factor of 'Equations 1*4 and - 15 are normalised error signals 
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and the third factor, containing the vectors a., and, b, are 
similar to the input vector x of Equation 4. Note however 
that this is only a, basic attempt to explain the nature of 
the vectors s and u and not a completely valid comparison - 
5 in particulalf~\^ vectors 

a and b ^represent the;^ backward ^ predictor 

coef f icienftsuof thfe^ElCT^ - arid are -Computed by 

means of- therrsample-CQva-riance^ matrix* R- by : * * - 

ap(f)=R^(t-lj 



and 







(16) 



.(17) 



10 jme ir^^rrejponding-. predictor;. eJ^or"^>wers can then be 

def in^aTu.aS:i. : i£: L : lt*j: „ ? „£E .r'H'-rZ fv:". i'{2iIJ 1 .;.7 . ... 



-~ 4^)=x°W^-[.xl(t>. 



and 



: u i 



c£ (0 =^(t-~P) -x: (t)- - -x^('f^^ ry •] 



- (18) 
- (19) 



1. As.Jl^J^Plicj^frpm above, the SFNTF 

operates- two™ "Separate- predict br—branehes- and by -enabling or 
15 disabling those branches we can create three different 
versions of .the ^algorithm: 



Version 1: Using f orw^rd ^ jmd ^j^kw^ra .predictors:, (not 
re c omme nde d ) 



Available values at tittie tY 
'20 From SFAEST 'algorithm 



.ime t : g//(< - if-and -y Nt p{t - 1) 



(20) 

(21) 
(22) 
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Verson 1 2: Using only forward predictors 

Available at time t: qjv(t - 1), jjv(i — 1) 
• From SFAEST algorithm: spi x (0/'4Wv r sp^i(t o ), 4(t 0 ),gp(i°),7p(* 0 ) 

- ^ : (23) 

S »M=[ jMtQ (25) 

9»(t) = 9N(t-l) + s 1 P+l (t)e f p (t)- S ^ l (t°)e f p (t 0 ) ' (26) 
IN A*) = 7p(«°) + <7.v(«) (2-) 

Version 3: Using only backward predictors 

•5 Available at time t: q^t - l), 5jV (t - l) 

From SFAEST algorithm: u F+1 (t), e P {t), u/> +1 (*°). e P {t°.), sp(*°),7H*°) 

[r]=-^-[^)] <*> 

g < v(t) = [ Z-l ] - w ™ 

J.v(t) = S :v(f - 1) r u^[(t)ej,(<) - u^}(J 0 )e^(i 0 ) (31) 
7/V.H*) = ip{t) ■+■ &v(<) 

Time update of the filter h v (t) : 

Available at time t: h. v (t - l),x,v(t - 1) 
New information: *(*).»(<) 



e ; vW = y(t)-^(l-l)x ( v(t-l) (32) 
,„(,) = J*fiL (33) 

h,v(t) = hiv(t - 1) - €iv(*)g*(0 (34) 

10 Predictor version 1 is more likely to diverge due to 

the fact that the calculation of 7 is not realised as a 
finite sum. Versions 2 and 3 include only a finite sum, of 
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P + 1 values in the computation of y . 

...As the algorithm of .Table 5 is rather complex and 
difficult to overview, we also include Table 6 which 
displays the variable occurences in the different equations 
5 of Table 5 and states which- variables are required for 
updating other variables.:. # ■ " •-•*."•' 



x Note: The definitions of 7 differ between FAEST and FNTF! 
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Table 5: Tasks of SFNTF> wdith' SFAEST in order of 



computation and number of required MUL/DIV operations 



SFNTF jwith SFAEST J , [ MUL j DIV 


Available at time t: 

-p{t - 1), a P (* - 1), h P {t r l)\ h jVf< -i, xp(t - 1>, 
7P(*-1), a£(*-l),a 6 P (*-l) 1 






New Information: 






Computation of the residuals and corrections: 

A.nriori forward / backward errors* 1 /residuals 1 — *~ •• - - - — 

1) eUt) =x(t)-&p(t-lj H x P (t-l) ; ' ' : 

2) e l P (t) = *(t - N) - b P (t -.-l)^xp.(i) ... _. ^ 1. 1 .. 


- - 

v P 




- 


Normalisation parameters • 

» ^'''-..(hm^^ ^- -, - 

4) 9 f (0 = ljL 7p+lW 4(i)(s^(l-l) + j4^i P (l - 1 )' > ) -. . 

aj fPW - an<1 /PI'; -.a-s^(t)x P (tj . i' 
6) Jb P (t - l) i \- p -f P (t - l)a£(* - 1) 


4 * 

i ■ 3.:. 

"-2 


1 
1 

1 


Difference of errors: , : ■ 
7) ^(0 = e^t) + M*-l)e^t) + AaV(«-l)g^(*-l); ' I < 






Difference of errors of correczed' filters" ' "*"" ------ - - - • --■ 

8) lp(t) - 1 ^ < , (1 _- rf((t)) ^ ( , fci ( t _ U ( 1 i.^p( t _ 1)) - 




1 


Corrected a-priori forward/ backward errors. ; . 

9) e P {t) = 4(t)-(i-7P!(t-i))^(«-i)p?p(t) : :• -•: 

10) aUt) =ja£(t - 1) + ^{f-m^y ' 7 "v ~ *' t~ — 

11) t> P {t) = l*,(t) - (i -TP(t))p&.(i) ' 1 • ; - : 

12) a* P (t) = - 1) + TPfO^m) 2 - - '— -• - 


~ 3^ 

3'; 


- 
- 


Time update 
Extended Kalrr 

13) SP+i.t = 


of the dual Kalman ^aih gp(t): 

lan gain vector: . . . 1„ , 

\ \ 0 1 'ci(t) ' [ : 1 ' 1 
'^p(t'r l)\ Aai(t^t) jva P (t - 1) ;•':<;' 


; p^i 


1 


Forward filter;;- ' 3 '■* - - — - • ■ •- — ... _ 

14) a F (*) = & P {t - i) -i*, P (t -xx/pinimt - i)p*p(t))&p(t-+-i) . 


""7"* ~ 




Dual Kalman.° 

w) r 


;a1n:~ f. ;. -J. 1 




- 


Backward filters = - i ' ' - ; * * ■ ; . ; 
16) bp(t) = bpXt:-l)± 1 7p(t)le% i ^ : 


; p ; i 3 




Computation jpf.dU^lXalman gain for SjFNTIF:. : j : ' 

17) see description of Version-.!)., 2) and Jl)._.: _ . ± ; . 


V 2 




Time update of the SFAEST filter h P (t): (not required for. SFNTF) 

ISa) ep(0^ Jt(t)-hf(feMJ. _ . _ ' ■ ' ' .■ 

19a) € P (t)±jp{i)ep(t) ' 7? : "T : ~ " : " 
20a) hp(t)fchp(t-l)-ft»&)gi>(t) , • 1 ; - : 


" P 
• 1 

-: 'P 




Time update of the SFNTF: filt^r-hA^v-' - _ . ' .. . ... ,„ 

18b) e^(t)> - h^. " - 1 ' 7 1 " : * ~ " I " i T ' 

19b) €i^(t);4= CAr(t)/7i3t) n > .; ' j i , : : > •; ' * 

20b) hyv. t = h^.t-i - €N{t)-Qtf~t ' ~ . 


■ r :iV " 

| ~' iv 


1 


Total number of multiplications and divisions: 


8P + 2iV + 39 
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APPENDIX IT 
Multiuser detectors for CDMA 

THE MMSE RECEIVER , ^ 

...Consider; the direct sequence spread spectrum" system shown 

5.; in Fig.-, 6,. _ . . - . ~v . *; . •:. 



* -In rit^we .haver M- tin = li.to M): users , each transmitting a 
~~ ^PN~co:de— oT T l~e^ ~"We~a^ume~that al~l~ the " 

users are : bpth :data' bit cand; chip ; synchronous and. that the data 
. -modulation: and the code ar take the values .±1 . The . signal 
10 at the input to the receiver at chip n is therefore':- 

y(n) = £ D m C WJ + W(n) 



{1} 



is the noise sample at chip n. Wiener filter theory 
[1] , [15] states that the optimal set of weights for an FIR 
filter h oot is given by:- 



Kpt = 0^0 yx 



{2} 



. 0yy is the autocorrelation matrix of the input signal and 
15 for na -stationary input: n : : • - =- . ' ■ : 



0 = 

— yy - 



:y{n)y(n - 1) rJ ' /(n). : . . -y(n)y(n> D 
y{n)yin r 2) : . _y(n)y{n-l) >;(/i) X (n - N + 2) 



. y(«>y<« - y(«)y<« -// + !)■ 



{3} 



where denotes the expected value of ~ We are interested 

in the FIR filter when it -contains no data transitions, ie 
when it only has one data .bit t from each user in it. Assuming 
lv that the data modulation* on • %a'ch code is independent, ie 
20 E (D a D h ) =0 , A*B by direct ^substitution of equation {l} into 
equation {3} or analogy with" the spatial radar case [16] , [17] 
, i ; rit^can fee . ghgwn ct:h.9;t * :. - ? ; ' .v. 7 

:.- ? -0^ = QAQ T + & H - : , {4} 

- :': r> rThe .matrix £)7has; .-dimension A&cM:-arid' -has * the- codes- as : its 
columns, ierr: r '-. ^;..:,:c .7 r. :. ; .-- \ t- . ;> 
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. ■ j2 r * 'denotes .*the...:tratnspose^- of matrix ; The - matrix A has 
dimensions MxM and has P m at position m along its "leading 
diagonal and is zero elsewhere. P m is the power at the 
receiver of user jtu. For : the ■ rest off ' ? -the section we shall 
5 assume P„ .-.= . 1 f or all:*m; -de perfect' po&.&r -c&atrol- with the 
power : of.- . each uaer normalised- :to "i r' The crM tertn 'is -the noise 
term assuming - the .-noise- - powers d-s- and- ; the noise is 

uncprrelated. \ : ^ - ^ ■ 



The vector 0^. is given by: 









yCn-i) 











{5) 



10 x(n) is the desired response vqwe S^e r .^nly -interested in 

the response v^hen the FIR filter contains the chips from one 
data bit, ie;when> n : = 0 and- at that point the desired response 
* is the data bit "for the" desired code .} Assuming the desired 
code is m =' 1 then x(n) = D 2 . Substituting this and equation 
15 {1} into equation {5} and assuming the the data symbols are 
•v , : uneorrelated. ^i/eiLds .• ri- '■ v.-.-- • . . '- .' ". -• 



* IN 



{6} 



Thus the vector 0 xy is the desired -:cOde . i-Pro'm -equations 
i 1 ) * {3} and {6}' and ; .under., the assumptions that we have 
stated, the optimum mean squared error performance is only 
•20.' dependent: on ;the -code., siet chosen:"" -The mini'mum imean squared 
error at the FIR filter output is given byi-v' . ■■=;.;/ ^ . : 
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..We Q an define the output^ signal to noise^ ratio for the 
filter as : - */ " 

Assuming that the channel bandwidth in Hz is equal to the 
chip rate in chips per second then we can, us.e the energy per 
. 5 — Bi t E h divided" by the " noise Tpqwer 7spec t ral~ densT ty N 0 \ as a 
processing gain indepenc|ei;t measure of the signal \ to noise 
ratio at the input to " our FIR filter . In pur _ case T- 



If* ~ 



Figs. 7, 8, 9 and 10 show the theoretical and simulated 

performance of the Wiener filter calculated as,aboye..f or a set 
10 of 7 and_ 31 chip Gold co4es 7 , v £ig.. 7 shows the, theoretical and 
_ fimulated perf oimance of .the 7 .chip Gold ".codes" in. terms of the 

output signal to ^o^se ratio from the Wiener filter. Fig. 7 . 

shows that the performance of ; the , optical filter with 4 users 
. . is practically, the ^me as the perf 'ormance u of a. matched filter 
is with no ,MAI Looking at the data, at ^10 dB ^output signal to 

noise ratio/ the difference between _ the.. $\ user curve, ^nd the 

matched fi.lter curve is ar purid „ 0 3 . dB . From. Fig.. .7.., .to, have 
t . 7 users in a systern . with -a processing, gain of 7, we mugt. have 

an increase of around l.8,_dB. in E h /N 0 . to. achieve. the same 
2 0 performance* as a matched^ filter with no MAI^.~ These- f i r gures are 

again measured ;,witb W^dB as r the ^ required' output .signal to 

noise ratio. " 



7 , ... 3V/ Fig .^8^ shows theoretical and simulated perf prmance..qf the 
3 ?- chip , G ?i d .^odes. - t^rms^ of the .output., .signal jta. noise 
2 5 rat:^ users,, we require around 

. - a - . i^creas^^ in^ E b /N 0 oyer; .-£he_, matched filter with no 

* v - M^i .... This ^ is £et te* „ than the 4 user,,. 7 chip, case , .partly 
because ji/y. is greater than f . Ift/.3JL B> When, there are 3.1 .users in 
:c _ , a C^ ; ^Mc^a S ^ S lo m 4.. lS* q^i^f d .ipcxeas^. ,1^ E h /N 0 is 1MB, over 
^.P'.the^ would, appear., :thajtl, the 31 
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chip system enjoys some inherent advantage over the 7 chip 
system, * probably ' due to the ' greater degree of statistical 
orthogonality between the 31 chip codes when compared with the 
7 chip case. However, other, factors such as the choice of the 
5 desired user or the choice of the 1 code sets cannot be ruled 
out . . r 

Figs'. 9' and' 10' show the simulated bit errbr^fate (BER) 
'"performance of "the '"7 "and " 21' "chip codes . .'These graphs show 
similar "results to Figs V 7 "and" 8 in "terms of the' increase in 
10 E b /N 0 required to* accommodate more "users ." A good approximation 
to Figs . 9 and 10 can be .derived from Figs. 7 and 8 assuming 
that because of the central"] limit theorem the noise output 
from the Wier.er filters is Gaussian. 

Adaptive "firtwf "receivers V '" " ' "*\ 
15 ~ In" mosf n cases " the DS-CDMA channel isTnon- stationary -arid 

as we'have" alre^d*/ 'stared "the MAI "will* vary "with" the" birth and 
death" "'of *^s icffta' Is . ^'Therefore' the FIR filter" used in -our 
receiver has tb be a' time varying^apprbximatipri to* the Wiener 
filter balcuia%ed"" using - an" adaptive"" algorithm We will 
20 consider- "the 'properties^ of the " twb most popular adaptive 
aigoritltifns^ ie"as't 4 ^mean squares '(LMSl and recursive " least 
squares' HRLST 11] ."We" shall * use the "1 chip Gold" code case as 
'" oU3r :; 'example , J aS'^the'^i ' chip xase" would* involve* considerable 
- eiktra"d6mpufeat?ion; i Fig.""* 11" shows the" convergence properties 
2% "b'f th-e^LMS^aig'cS^ithffi with' ttfb different values' of the' adaption 
' Ratafee tet* "/I. * Thi^'gr'aplis show" "the feedback - error ensemble 
summed over 100 independent trials. With '/T~'=~ 0".007 the 
algorithm converges to 20% of its^ initial value within 
"a'pprok'iitiatSly ' 20' "data bits , * l>uif "wit'li 1 "ju "= ''6 70007" convergence 
30 takes" acsrbund" i : t)0 J "dicta b"its"V'"Tri a ^ "DS'-"CI)MA ~ cellular" system, 
whe're^the ^charihfel** varies * fela^ively/ quickly' with" respect" to 
^ the Sat aerate b*f ' l & 4 sf ngfle'tiser/** this' slow' convergence ' will not 
~"bfe ab'creptabie . v Figf ~I2 shbws t*Ke "BEST perf orimarice "of the LMS 
J "algbritfim n after- coiitrergWhc% with' the' same values ~oif fl as for 
•35 F~ig; r l"r. It shows" niii&r - the" "smaller valufe' o'f :i ^produces a good 
' c approximation t o" t h & Wien^f'f ilterV*6ut the perf ormance" of the 
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LMS algorithm with the larger value of /i is around 2 dB worse 
in terms of the required signal to xibise ratio for a* bit error 
rate of 10**. Thus it is difficult to find a value of > which 
will provide a fast "enough convergence without introducing a 
5 large residual,, eVr or into the performance of the filter, even 
jor the^7 chip case . The convergence time' wilT'be 'greater if 
the codes are longer, if we consider a DS-CDMA cellular system 

with imperfect power control and a multip^th'channel , then_the : 

~ ~ eigenvalue spread of the input signal will be, increased. This 

10 will also have an adverse effect on the convergence time of 
the.LMS algorithm. More evidence of this is contained "in [18] . 
Thus, the LMS . algorithm Is" unsuitable for this application. 
Also shown in Fig' li" is the convergence of an RLS" algorithm 

. . /with X = 1 and the "diagonal elements of the autocorrelation 

15 matrix initially set "to 0 . 001 . ' This shows typical behaviour 
for an RLS algorithm, rapid convergence by The : RLS curve 

in Fig. 12 shows ^that the RLS ^ algorithm does not add any 
significant residual error, to the Wiener filter solution. Thus 
the RLS algorithm^ is potentially more* suited' to ~" a DSCDMA 

2 0. cellular system.^ _ , . ; . . ' ^ 

Discussion ~ . „ _ ~ ^ _ * 

In this section, we shall look at^ the assumptions made in 
the derivation of the filters and look at applying this type 
"of receiver to .a cellular "system. ^ ' _ ~ _ 

25 We^have as turned that, all the channels are data bit and 

chip synchronous.. I f Jth^ channels ^re not. chip synchronous, 
.. provided the . t receiver /"samples _the desired . channel 
. synchronously, any^ change^ is likely _ to be benef iqiai as the 
effective.^ the^ channers are "not data 

[ 30 bit ^synchronous however; „the effect is^ greatly detrimental. 
With transitions pf^ fcSe_ interfering channels occurring in the 
^.middle of the desired^ channei data bits , each interfering 
, channel _re.quires two eigenvectors in, the ^ Q jfiatrix. Thus 
performance, is reduced. and the ^system breaks 'down when the 
v 35 number of users exceeds 0.5N instead of N for' tjie data bit 
synchronous case . ~An interesting case is when the data bits 
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are. almost synchronous, for "example, as .synchronous as the 
transmission delays in the system will 7 ' allow. In - this case, 
it may be possible to ensure that any transition occurs a 
small r number of chips (compared with N) from the beginning or 
5 end of the^ desired user's data /" bit.' In "this" case" the 
detrimental effect may "not "be ''"so great, although' "this 
hypothesis remains^ u n P^°y en - We have alsb " assumed perfect 
power control . There is ' some' evidence^ that Chris type of 
adaptive filter structure is * tolerant to variations in the 

10 power of /the received codes"! 19J\" The ^last assumption we' made 
is that the data' bits are uncprreiated.'" Provided k 'tRat the 
signals, are , independently generated/ '"idle"' channels," are 
suppressed and that care is taken over the content arid" timing 
of. "control information/* ^ this ' should /^be ' a~ reasonable 

15 assumption. \ , , ~ ~' ^ " 



If we are to apply this receiver" structure to* a cellular 
BS-CDMA system "we heed" to take ihtb ac^courit the' ndn- stationary 
multipath channel " and the birth and" deatfi* of "users . To take 
into account the multipath channel, we can lisre thie" technique 

20 in [20] , replacing the impulse response of the spreading 
process with B^, where is the convolution of the 

" impulse Response* of "* the " channel/ with ' the "impulse response of 
the* spreadiiig process and repeat /tHe "analysis. 'This paper 
already describes the theoretical optimal with multipath 

25 evaluated for all users of the system simultaneously. The 
non- stationary elements^ bf~*tfieT"cha^ taken into 

accouht by /blocks ~ 6f training data ~wit£'~6r~ without '-decision 
f eedback in between, /MeT birth" and' death ^of signals problem 
will * be - greatly alleviated" Iff ^ll^usfers are ' constrained to 

3 0 switch on ' and' ' ofl'"only ' lmmedia'tery proceeding" a " training 
block. The Convergence 

algorithm and speed'at"*^ 'determine 
tTfeTlengt^ £locKs . Hovfever, by 

employing* the "RLiS af gorithm^(or ttie cbvVrxance"f brm"6f the 
35 "least squares "al go rit Km on training*' blocks" [ IT )'/ tlie ratio of 
training data to Information carrying 'data' can '"'be kept 
reasonable/ * " " ** : -v " L ' " — * 
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-Sbme" desirable properties ^cf this *:: receiver are: that it 
• '- does riot require to Jcriow" the --desired or ' the interferers 

- spreading codes, provided "the .training data is known and its 
adaptive nature will- allow' it- to ^reduce the effect of strong 

5 -interferers" from neighbouring" cells and local narrow band 
interference-. - : 

- Conclusions - 'il\:\^~\ : 2 J\ : '~\_ : '''' 'il '■* S" 1 __ __ " 

We have shown that an FIR filter can be used. to separate 

a desired signal from MAI in a DS-CDMA system with only a 

10 small degradation^i-n Gaussian noise' performance . -These filters 

. w :.-can" be approximated; by- -trained adaptive -filters. These 

- sub- optimal faltei-S^may make practical receivers "for- a DS-CDMA 
-~ cellular- sys t err. . -r.:; .. r r- V. r .c 

- THE DEGORRElATiNG - RECEIVER : * ; -~s 

0.5- *~ : r : -The^ de-dor* ei^ting ' ^d£te~cfeor - is - similar * -to" --the- MMSE 
■ 7 -"' detector except- chatty the 'no is e c tertn J is-Tiat taken "into account 
ie. the- : s£coi*d -t^rto -ili'eqa^tioii -{4} " abov^" i-s* he^l^ected- in the 
- ^ f orfrtul aS s f or' ^ &ffr - ~ ?■« ' - ■** - r - - ' • • - ' 

- - "RADIAL BASF'S- FUNCTION -RECEIVERS ; v~ - - - ~ ^ Z - 

20 :: - ' Introduction; : In - many -DS-GD J MA - communications systems 
there "are three" ?sbur-££s of distortion -when^the -signals-arrives 
at the -receiver >" j Struetured-multip^ inte'rf ereiic : e : (MAI) 

froth other users -in' the-'syst^m' r Gaussian- noise -which can often 
be extended to include unstructured interference ' knd time 

23 dispersion due to. mul^tipap^i ; (prppagation. The simplest DS-CDMA 
receiver structures are based on matched filters for a 

' * noh- dispersive r dha^ei-'^ -multipath 

- • c^ahnel^V Th§ -^'^ Wtched" fil ter /RAKE ^receiver s 
: ' - can- -be" • ehnanced- by applying canceiratioii" at ' the expense of 

30 increased receiver crompXexit-y- [12] . Many- proposed" receiver 
structures are based on linear equaliser structures. Examples 
include- dec'Q'rrelating-- 1 " receivers on the 

y z^ro-'forcirig ".eSjaali-Ser"- - and" those- based" on -the MMSE equaliser 

- ~[2iK*V ; W^ receiver -stnictures'with the 

- "3 5' ~RBF :rr -receiver- . " Thi-s :: *- nGril inear - receiver minimises the 
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probability x>£- ecro^ when^ deciding ,Qn a ■ data -..bit . It has 
already .been . shown: - tiiat . -an . RBF equaliser has superior 
performance to a linear,* equaliser in multipath. channels for 
conventional minimum bandwidth signal ling/, [2 2-3 and that -an RBF 
5 filter; Has. -superior. : performance -to: a --MMSE : filter , % JLn a 
non-dispersive CDMA system [23] . We shall show that the RBF 
based receivers have the ability to considerably increase the 
capacity of a DS - CDMA system when compared with -other receiver 
structures., - ' \ ^ . , r m ~ C} 

10 * - AWGN jchanneX sy,steih. iiiadgl:;,: -To. enable.- us- tp v compare -the 
•■" wide variety of; .receives, .discus seA . above y,e' shall first 
... • consider -.aix^additive ^white: Gaussian -noise .Chanel ; .Our . system 
will consist of U independent users each;- transmitting a 
DS-CDMA signal which is chip and bit synchronous, and with .all 
15 users transmitting equal power;:.no.rmaiisedT.tQ:.i :i;Th^ data .bit 
transmitted ibyrjuser; u du-ring?rbit^ be ,denoted>by 

P u ,(^L:;and tfee spreading- code -£ o:io^e^ by >q^^ We .-will; use n 
to de^evthei cMp^nu^^r-withi^ tfce .cp# which ~is_an , integer 
between 0 and {N~D where N is the spreading*:. sequence -length 
20 (processing gain) . The spreading codes used throughout are 7 
chips long and randomly gener.$t;eta . r Wrt;fto.ut; loss,- of. generality 
= - . we eanj-assurae^that ;the desire^ uger : is ; >user ,0 r . .The channel 
:.;mpdeA^ ^white ^Gaussian noise. (AWGN) 

,*. m^dei G(kN + i) . Thus 

25;. t-hft; chip v rate; sig^ -denoted Y(kN 

z "tJiA ) iz ifilri ibe.v- '* -thv ii;:: 1 z:j r. ?. - . r ; 

- ~ .j^-i^isr* r:~> r-s^^a szs 3 -3 y.:.;z -jir: ... 2 z^- - \. 

- . at-thecEoint inhere, received* -In r the AWGN 
^ase- there^,isj : n^ range w «, 0zi<N 

^^foxoutside 3 this-^i™e n£ the ? . rS^gx^l ^ili j:r c<^tai%ano c -useful 
3j0 o inf orgnat^oji relating :to dat^iifeote. Tev.^:^': f. 5r v: - - : - 
t-t,:<Z • r-:;* r.„r- z Zs.: *:.~^rz~ :zc rs^i.-rl e^c. ■ •? 

— La&FGN c^ajxqel : receiyeir : s^uct5i^^ ; .Vter ^h^ll consider three 
Jti^j^:peceiver ^ rthe^MMSE receiver 

-anci. ttte- -.RBF- «oeiverr lntall r >^s^&qw^ c shali 2 ajs^ume' -that the 
,.cod§:S:.and the .§dgna : l powers l^re known; at .^.^ceiver ^ -In fthe 
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. ..non-dispersive case the ; matched: f ilter receiver i's . given by 
....an N.. tap. FIR filter whose .-c or efficients are the. spreading 
. code, for the desired receiver -_ - ■-_ 

The MMSE, receiver - is also an N . tap. FIR filter, but the 
3 -co-ef if icients of this filter are ..given ( in'v vector :; form) by 

' [24] ; - 7 i. • ./. • : ; , - ■■; : -.. .- .. 

~" V~: _ ." ~ ~ M™ 5 ^ ^ ; - -• •• : - •• • • ~"~T 

is the autocorrelation matrix of the cyclostationary input 
signal with dimensions NxN. 4^ i s the cross correlation 
: -••:. ve . c, r: <:5r • :.7 h ? receiver, is. a nonlinear filter whose estimate 
0 of ..the- data, output la^gi-ven-by.: rr. : . \:.\. . 



Ml 



where. ^(Jc)_- denotes -the-; input: vector, -consisting of the-N chip 
„sptaced input samples ...at .=,;dat a,;; bit;, time-" - -*- ; and- c; are the 
centres. The centres are, -the noise- free' input vectors f6r all 
possible input data bit combinations. There are therefore n. 
= 2- centres, j . | denotes the 'length' of the enclosed vector 
(Euclidean norm) . a is the standard deviation of vthe, noise . 
w. is the value of* D 0 associated with centre "c ; . In our 
simulations we shall assume that the centres and weights are 
Jcnown.at the receiver-:,; For -all ir.hree receiver' cases we 'shall 
a l s P consider^ supplementing :the -.receiver: with a-J : s ingle '-stage 
of -parallel cancellation -similar' to- i[12i]-\ -' ■ ■ ' 



o -1 



AWGN channel simula.t^pxr:r,esul^s>':The"xec'eiver -structures 
above were simulated using Monte-Carlo simulation. The graphs 
:■: - sn ?. w :^. e ^PSii PER- -averaged ^0\*er-^ all ^- active users 1 in 

r . t ^ e r^^? e ™:-.P 1 ^- t . te . d - against the- -.number of active users. Fig. 
1.3 shows. resul.ts, : f or E^/Ng r_=, 9: dB'. ' These figures .show ^thait the 
jR e fF: f ? rTO - anc f : of J. m .at?hed filtering' be-cbmes"" : -poo^ as the- MA! (the 
. jiumber , of ra-ctAve. £ users>..:.increase&. .'The- nuinber of - act ive : users 
- :¥$^? h ~&yL BS .-an .acceptable: BER* for..:a':-matohed ; f ilter receiver 
?.: considerably.- .*tp&.. the -.adtai^iori- "of - cancellation , as 
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.this filter :is clearly interference limited: -T'he MMSE receiver 
.. performs better than: both' the matched '•'filter' arid' the matched 
filter with cancellation, -.it too improves with the -addition 
of cancellation, although the - improvement is not as great as 
5 adding cancellation to a matched filter. The RBF receiver is 
considerably -better .than both^the matched' filter and the MMSE 
receiver,- . with or hwithout- i cancellation. However,- the RBF 
- receiver is not improved by the addition of cancellation. This 
is because in this scenario^. the/ RBF'- is the maximum likelihood 
10 receiver and therefore cancellation cannot improve its 
-..performance:.. • -.v „--':; -. .: - : r . • 

: . •• , * Mul tipath -channel . .system model / -We^ shal 2 : now c'o'n^ider a 
stationary multipath channel " model . '■' The 'channel '* we' 7 "shall 
consider has impulse response :- 
15 M(z) = 0.3482- + 0:'8704"2^ > '0--3482Z- 2 

All signals are assumed to pass through the same channel as 
would be th£.- c&s.e: -in the . downi3.n]c' df a %iobile r ad fo : system. 
The r^ecejiyed signal at. the titoe ^ata^b£t- : 3c, ch'% i is "received 
by the i direct; path and therefore^ -becomes 1 ': - ; • " : "* 

■ * • : - ~ . • ' (u*\~ ■'• - - " \ '-' '-'^ v.: «.;••.:..• - 

Y(.kN + /) = 0. 3482 £ D^Cuj + P(kfi +. i) + . _ . . ; 

• &"■ ':. 3 -.- ' " • - \*«-o* - ■- - -...-j - -•---/ '-" 

6. tfU^buWC^^Gdw Si - 1)J + 0.*3'482f^; 1 O^kyC^ + C(kN + i -2)1 

. .2-9^ -NQAe that; i& i-x..ds negative '-the-n- the J> u (*)-£„.£•, is "replaced by 
- . Land; J.-fc.-j <o.r ri.-.?c lis- greater - : -t^air -(if - we ' : extend 

the receiver .filter, v length.: vto&yom^ -2*)'* -then' i£(*Vc u ;r.* is 
replaced by Z3 L . (k+1) C u>i . N . x . This is a direct illustration of ISI 
, . caused -by- the t multipaith; ■channels. '. j'-i -:. "„•&»:.: -o-.:. : 'Z" ' 
-.:.,.rv. r-T ; - 3 ■-• 'tvK^J'r: . z ■ ■■.^£- 

25. Wu^-tapaffc .. channel, irece-iver?- structure's t . We-' shall base all 

ou F. , n .f w receiver structures ^n'-arhevF^et-^of '• input ''"'"signal 
r: s - ai ^P le : s -yiWt i )• - whexe" the ; range - ©& i" I ; s -'ext"4fided' to"" oii<W+2 
.rf. ^.capture %il,j:he signal? ehe-rgy~that''origanated r from data bit 
-'v j^-r^his.- combined with the ■•■mult±p , ath'-Ghannel%"£4l\-esu=lt : -in ISI 
?°- f^P" 1 both • the ^recedijig. ahd f ollowihg-da^ta bit ?-Tlfe equivalent 
i .•> P^.'r c fef : ma P.?& e< 4 -.filter $®e a.. muli;ipath ^haiHiel^i's aif F^I k" ; "f il ter 
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with N+2 taps where the taps are given by:- 

The * represents convolution"' and' : .c 0 and M represent the 
spreading code of user 0 and impulse response of the channel 
respectively. The MMSE receiver coefficients "are given by the 
■ 5 same -equation'^ as- before/ 'however, 'the "elements of' the 
( (N+2 : )ii(N+2) ) -autocorrelation matrix ^" and" the" (n+2)~ cross 
-.- 77.:— ccl:rei _ ati^ri'"'-Ve"ct'or"" tm-se^be - ^alc%l : ali^"d "acceding*" "to' the 

• " '" derivation ' of -thS- v &MSE- received ~ taking info " account the 
multipathr- The' -rbf-"'-- receiver " : equation 'is" also" essentially 

: 10 unchanged/ -"except ' now'"' the input' and- centre vectors'' "are' of 
'- " - length ~N+ 2- arid -the ::r number - bf centres : " il " 'now?"' 2 il ' as 'all 
possible combinations of "the previous; curre'nt' and" next data 
bit must be considered. This number of centres increases 
rapidly with the number""" Mf ' 'Active ^Users' "^t/"' and the 
r 1-5- 'computational'" load will rapidly" become " impractical. 

"" \ ; • Wul-tipVth ctianhei s'imula ti on 1 res62 ts~:~ Gr ap"hs' Tor' the BER 

performance bf'^the" receiver structures" " described ? in the 

- ■■ - -'previous- Wctiofi-^re 3 shown* in • Fig 1 : ''x4 for ~E b /N 0 ~ = 9 -'dB? These 
-■" 4 ; graphs'- show- -simlar trends to 1 th"e : nbrtdisper^ive case". The RAKE 
' 20 receiver' rapid! ^"breaks ''down a's the MAT interference increases 
because the interfering users" are r n"6t tikWn "account of " at all 

• •-" - in the- RAKE receiver "and are'VheVef ore 1 ' treated as" Unstructured 
-'- ••- noise: The MMSE i "receiver '" does ; considerably 1 " be^teV 'than the 

RAKE receiver : The " RBF- receiver pWf or^' ^etftif than either 
25 the RAKE receiver or the MMSE receiver. However, it does 
-" - - involve 'a considerable 'ifecYease "'in computational "complexity . 
the ■ r RBF ^"f esults 0 '- ari c "t-runcitfed" ' ; a% s "'^ ~ Users'* because the 
-•'cai-cul-ataon-of ; '-e : 'users' woufd -require "eVaTluating the Euclidean 
" ' distance f rom J 2 £8 *' centred for' every' data' bit sent."" '" * ' 

-• 30 ~ a Sisbasiion- and-Conci'iisionsV We have "shown tha't ah RBF 
- -" receiver s has -greatly "-improved "perf dr\nance over . the more 
----- '• conventional-' 2 matched"'- f-il : te"f "and- " MMSE" ' based ' receiver 

• ^ ■ • -structures. Thfs performance' increase is apparent in "both AWGN 

- and ^ muli"ipiiBx i: ch ? annels 1 . ;? However "5he" computational 7 complexity 



BNSDOCIO: <WO 0051260A1_L> 



WO 00/51260 PCT/GBOO/00649 

- 37 - Appendix II 

and the number of variable parameters in a non- stationary 
environment presently make- the RBF filter receiver impractical 
for mobile radio app*Licatipns . 

NEAR OPTIMUM RECEIVERS 

5 . Neaf -Optimum receivers use similar methods, to the Radial 
Basis Function receiver:, usually .with, simplifications to the 
Radial-. Basis Func.t i on- „by dropping the. exponential, term and 
simplifying the. _ Eucjlidea"iX-. r dif £ erenpe - expression. These 
receiver^ - can . be. _ implemented, ,chip» level - or the bit 
10. .level _ (af ter 4 preprocessing) _and . /often.,, involve . the Viterbi 
forward programming algorithm to search £or the- optimum centre 
(path in^ t^e viterbi algorithm)^ ... . . . . ..... . ; ; ^ 

CANCELLATION BASED .RECEIVERS ,- T . . - v -u ■. r 

The^optimal ^recQiyer- f ar, a di re c£ £ sequence ..code.^divi-sion 
15 multiple access (DS-CDMA) system is the maximum likelihood 
sequence, r estijaatp|:. _ [25] Jl4] r .whiqh^i-a -ef^ct-iyely an RBF 
receiver witl; ^Jixif inite rnemory^.. Hpwever, ^this receiver's 
. complexity rises r exponent ially D with... the-^numbjei^ of - users and 
_is therefore .impractica.l-^ At the. other- end., qf-^the complexity 
20 scale, _ a matched .filter .is ^commonly... used:., as-. the receiver in 
a DS-CDMA .system. 4 A ma^che.d. .filter is the optimal . receiver 
_ only in additive whi^e Ga^s..sian #Qise^.(AWGN.), and has very poor 
performance when -tlie. lev%l^of.. ^mu^t iple : access interference 
(MAI) f ^om^pther users iji. th^e _ system, is high... .. - - - 

25 Between _ these two. o ^xp^epesr a, ;r range., of . sub:-, optimal 

receiver, structures, hayf^ fc^gn^pr opposed , The^ma j 9?rity .pf these 
. .car* bje split. ^ are 
typicalJLy; based_ on mat phed..,3E altering, followed by r subtraction 
of the matched filter output from a delayed version of the 
30 input signal [12].. [?6J [27J .?c^a^ser.,s trxLptur.es-. use a linear 
or decision f eedback . FIR, ..filter - -with ? the ...co-efficients 
^ optimised according to va.^ing criterion su<:h.,as minimum mean 
. . square t errpr t (MMSE)_ _pr^ ^ zjsrp-f orcing-. [JLBXl20hL2Ji] ^ A 
comprehensive review : qf, this ,area -a^ici.^ ref.er.eixce. list, can be 
35 found in [21] . " 
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In this section the performance., of a Wiener, filter 
receiver is compared with that of a simple parallel canceller. 
The combination of the two- -techniques is also examined. 

System description -r 

5 The system to be considered consists of U independent 

equal power users transmitting , both data... bit and chip 
synchronously as showi} r ip : Fig.,. 15 . The data bit transmitted 
- b Y --us.er u_at- bi t. ^ime--^ will - be : - denoted -by - D^k)- and— the 
spreading code for user u.^will be denoted by, C u;a . ,n denotes 
10 the ; chip within the code , /tyhich is,an integer -between . 0 and 
(JST-l) and N is _ the spreading. ^ sequence, length (processing 
,.? ain) • T he spreading codes used an ..example, throughout < this 
j Pape^ are ^64 chips . long apd_ randomly : . generate^- ^Orthogonal 
^ w ^ s ^ P r P!?^rPr%hqgonaX .JGold) codes ., would give better 
15 perf ormance, T b^t would, require longer . simul at ipn,, times to give 
. mean ; ingf : ul BER results . . In .many communications systems the 
channel will destroy the or thogonality -.property , of orthogonal 
codes anyway. Throughout: 'this' paper the chip and data values 
will be normalised ' tD^ ; "^Withoift 'loss "o^ generality we can 
2 0 assume that the de.sir.e4. .user is . user- 0-^ The channel -under 
consideration will be ^ simple. AWGN : clj%nnel, ? with, :: the ; .-noise 
;^ eri oted by G{kN+n)_ having., variance a.%,.,Thus the chip.,, rate 
sis™* 1 arriving at., th^ .pceiyer, ^denoted-* YlkN+n> :J • will r be : - 



. Receiver structures and .theoretical performance - - - 

25 The four receiver structures, .shovm- . in r Fig . : - 16 rrjwill be 

examined, Receiver structure a) is the simple matched filter 

. . - ca ? e % Ti > e - rBce^yejr.-CQ^sist^ of : _an : N tap -KIR ; filter whose 

c ?7- ef ?f S^ en £f -A* ^e_^ :r sqaled version .of : the_ desired, user' s 
code , ie : r , \ ; , . _ . ; . ^ . ^ r ^ . . = ..... ^ . . . : . . 



30 The matched filter treats the MAI as noise and therefore, if 
the signal power of each user is one .and the codes are random 
the central limit theorem -allows us to approximate the BER for 
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- • (1) 



where erfc() is the complimentary error function. 

Receiver' structure bj " 'is "a 1 simple'^ "parallel canceller 
using matched filters ■ ta despread each interfering signal and 
5 reconstruct ■ a :/ replrca ~ which is "subtracted frotn an 
appropriately delayed : Input "signal, daricelling 'much of the 
interference , The remaining 'signaf i v s**then "despread using a 
" matched filter to the* desired ii£er. The reconstruction of the 
signal ; assumes that "the poweV of ^the -signal"' Is known exactly. 
10 -Thefefore^afi^Iriterfering' signal which is "is received" with the 

- correct sign will be 1 cancelled * exactly " and ah interfering 

- ™ signal " received "ihcdfrectly" will " have ' its 'amplitude doiibled 
* <powef : quadx^pied) : '->- A 'lowet ' bbUhd* on 'the BER £br tfiis 



structure -is* ^gi van - by •? : 



B£/? = erfc 



1-5 -This" equation is only * a ioweV 'bouhd " f 6t £wo reasons.* Firstly 
~- the- noise "samp l«s* -for- each matche'd filter 'in* the cancellation 
■- stage are the- -same, "Whereas the "above equation assumes they 
- --are- independent 1 - ^ampies Secondly , after "cancellation the 
combined interference plus noise distribution of the remaining 
20 signal is a poor appr6xirtiatioh£t6 Gaussian because it consists 
of the Gaussian noise and a limited number of enlarged 
interfer&Sfiig 1 : s i'grial s - which' ~ Waif 1 '-fife' a"*p6or*' approximation to 
-a Gaussian 'di3 : tribut i oh - L ~ - f ~-- 



• "-Receiver structure : c) ; iV* th^ ' Winner revinson "or MMSE 
25 - -producing " f i 16 e r - [ 1 5 : h frr] V It ~£s "also^lh W'ta^'FlR" filter',' but 
the co-efficients of this filter are given (in Vector ' form) 
by [24] :- 



Wh'ere- 



h = <&~ l <b 
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The T superscript denoted a matrix transpose. * w 
^autocorrelation matrix of ^ the' ^input -signal with "dimensions 
- -N^CN," given- by: ■» ! ~ . . 

<- -~ ■ -' - ®„ = QQ T - + * 2 I * - : • 

: : The matrix:* ^ *iias : . dimensions -JficD"- and has t; the ^cod'es as its 
' 5. ^columns-, 



x£ : - v 

^•0,0 

2 = Co - 2 



Q,0 



^2.0 

-.CVi 



■ r ■' Cu-12. • 



.C6m-\ Cz.n-1 • _ 

J is the JV&tW identity matrix. is the cross correlation 
: vector : -f : .:• v. 1 - : r .1 .-?.r? - •• ~ -r- ■: •: •_ * 

•:The MMSE'. is-'givem by:-" ~C v:.i u ■_-„■;-■-.: . :, .= 



and the ; " BER? performance ' "(def i vat-ion- A) is 



BER Wienfr = erfc 



10 




A^,)2 * (2,-0)h T (p yx - 1 . 0+..MMSE- 



15 



To. ;obtadn^^avera*g^: 'BER^ performance^ this- -ecjuati6h : must be 
-^averaged' tjve'r^&ll* ^tai^^if dear's"-' r i-h the- syst*em ; to* give"" BER^ e ^ . 

Receiver 'Structure ; d) - -xs- a par-a^'ei- ^caneeiieif- ■'■ using 
.Wiener miters- fo'r^the : irix^a^ 

data bits. By analogy with- equation -2->.--a : lower -bound;" on the 
BER performance of this structure is given by:- 



BER = erfc 



2 0' 



Results.", v.r.s- : . \ ..~r\.-~ *■ ?.-■- : ~r. < ■>*" 

: " --All- simulated' 1resuitkiin%his" sect ion'"' aire 1 f of "W/OCK) data 

bits , - .with sfeKe -BER- averdged i : over " : al£ : - usersv Graphs show the 

average *BER al- r l users in- the sys^m' plottecf against* 'number 
.6fi:u'sets:. v Th<§ .'sequence -.leiVgt-h : as -647 : The'~signai * tb : 'Gaussian 

nois^ *rati"cT is- expressed- as ' r :: - • .' - 



E b /N 0 = 



2a J 



FigV -1'7- Show's-- a- comparison- between-' the theoretical 
- .'result derived -itf -abcwe arid simulation results. The : * "results 
^'how- ^ fehafe-^ egtiatiohs- (l>" - ^and : (3^ -give -good 1 ' and' 'Unbiased 
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estimates of the simulated^results fLo.r the matched filter and 
Wiener filter respectively. The deviation of • -the : simulated 
results for the matched filter-- from equation (1) with a low 
number : -of users - is .caused by -the interference distribution 
5 only being approximately Gaussian ,,and t:he limited /number* of 
simulated points. Equations (2) and (4)- however only give an 
approximation to- . the .actual -^performance of the two 
cancellation receivers. This estimate is biased towards a 
lower BER for reasons that' were discussed above. 

10 Fig. 18 shows simulated results for a range- of r signal to 

noise ratios. These results show- that 3t matched filter is only 
a good receiver structure when Gaussicim=noise .is... the dominant 
source of interference. This^only ^applies when there are very 
few interfering users and/ the r - : signajL;; to * noise.- ratio .is^ low . 

; , ; t • ; \ -•- - ' '* 

15 A single : £s : tage"*^ on matched 

- f ilters . -performs _~ better.- than.,: a.-;.inatg:hed.v: /fiiltJSjrroexcBpt' ."at 
extremely high „levels. : of : inte^f^r^p.e / ; .>rtia^e';;the JJER-'is r . too 
low for most communications systems anyway. The improvement 
. is r fairly ?. s with v £^/^;.-?=v,.9- vdB,jathe; cancellation 

20, receiver ..cait support _ three jtimes as. many u^serrs; a-s the -matched 
. ..filter if ;.the required BER is. .ljG, :2 ; -y / : 

The Wiener fitter will- always perform better than the 
matched filter. The'"bnly)~c^ of these 

two receiver systems is the same is when there, .is:. =210 MAI 
25 -interference .and^-ti^ Wiener^ fdlt^^-Jp^jsepa^a^ra ^scaled version 
of r the, ... :ma t c^ieA. filter .^-jrtite.Y;. .Wier^er^rif ilfeer ;; .-also. ? performs 
significantly. . ■ better..^ than^ pite^i: matched.-;; filter-. -/-based 
cancellation ^ receiver . : yi±th : ^ b /LN 0 S^^B^and : jac- required- BER "of 
10" 2 , the Wiener filter will. : al low apparo^ as 
3 0 many users as a matched filter^ and twice as many as the 
matched filter based cancellation receiver. The computational 
.coir)p.lexity of : _-the r ,Wiener ; : f ilt.er receive^ is.: identical to the 
r . _ m^tchred ..filter, receiver,,, and ^i~-3^^^^^^^.E-2' e ^ s ^'^ s ^ tiie 
, ; .matched filter based cpcellation. reeeiyer- ; .r This^assumes that 
35 the Wiener co-efficients are calculated in advance. 
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The Wiener filter based cancellation receiver gives the 
best performance of all . It will allow approximately seven 
times the number of users if E h / N 0 = 9 dB and the' 'required BER 
is i "10"V compared with the matched 'filter receiver. 

5 Note that in all cases there i*s a significant variation 

in performance between users/ For example, using the Wiener 
_f Li'lter *^receiyer* : with 48 users and E h /N 0 = 9 "dB, - the average 
BER is 0 . 00"940' but the' user with the* best spreading' sequence 
expediences a BER of 0 . GQll&i * arid the ' user with the worst 
10 spreading / sequence - experiences ' a " BER "of 0.038092-. This 
variation is dtfe : to' 'the '"cross 'correlation properties of the 
code set and will also apply to an orthogonal * cod£ set' which 
has lost its orthogonality due to a'" multipath channel. 

Discussion and Conclusions' 

l5 "•' in this ' -paper J we laave" showii" that cancellation' and Wiener 
filtering both perfbrA"mueh better /jLnf MAI" than a simple DMF. 
Wiener fillers in particular of fer.-,a* large increase in the 
' "rxumber; oj&i i u^rs- ^tha'fc .-c : atrf b£ accommodated. The two ideas can 
be successfully combined and provide a better performance than 

2 0 either on its own. *" 

— - There is : a "cost to* pay for improved performance ." 'The 

cancellation' -receiver "introduces : a delay of one data'bit and 

- requires a substantial * iifcreaie' in either computation or 
hardware ; * Both *the* Wiener- filter and* th^ cancellation ' receiver 

25 require knowledge of the *flurti>er arid"" sp^eiading sequences of the 
interfering. users.-at^thC receiver, and : the -Wiener filter also 
requires a priori knowledge of the signal to noise ratio. If 
the signal ..tp ^noi^e, ratio vis. -very- high/, calculation of the 
inverse of the autocorrelation matrix for the Wiener filter 
''3 0 "can be" • difficult * 'as *it : becomes ' close to singular. In a 
compute'r 'network*" 'application the calculation of "the Wiener 
" ' " f iit6ar in" "advance * "may - r be feasible , ' but in a" * mobile 
* 'CommunicaLtibnW "appl'icatibn the" signal" to noise ratio and the 
\ number of users varies rapidly with time. However/ a good 
35 approximation to the Wiener filter can often be obtained using 
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an adaptive algorithm ^ 

Derivation A . . _ . 

T9,.calculate- the J3ER -for the Wiener filter we require to 
find a relationship between the MMSE and the variance at the 
5 output of, .the filter. 

. .MMSE. = E[(Z?o(« " «?1 = E[D%tk)] - (2. 0)ELD a (k)x] + Elf) . ^ 
whereE[] - denotes the : expected value. and ,x denotes the filter 
Qutput . ^D 0 Uf) t^es.. ..the values. ,±1 f therefore :. r - : : . 
. \ [ . _.\ MMSE = 1 - (i0)^ T 3>, x + Et^] " ' " V ^./^ JA.1] * 

The variance, of the filler ,pn,tput, is. .given by: - , : . 

where Jc denotes E±x] r „ _ . - 

10 Subtracting [A.i] from [A. 2] and rearranging gives : - 
var(i) = - {h T O yx ) z + (2. 0)h T ^^ 1.0 + .MMSE, . . ^ , 

T r? re f:- r 5^. < ^s.uming. .,th^. ._qu£p^p ; _ ..of the , Wiener filter has a 
.Gaus^ BE^ ; ..i^giv§n by;- . ; v ^ 



_, >.,.-., . .A . ^ r .; 




*C2..0)fr J, 45 r - 1:0 + "MMSE' 



VOLTERRA RECEIVER 

The Volterra receiver is made up by a power series 
15 * e *W*P^^° n ^ pf -i 1 ^^- ^peived . signal followed by a linear filter 
j.:^ u Sr .v? s - ** e ?5 th^, . Ypl^jrxa ..expanded signal v. is used 
. K . : in^stead, of . ti^e r origina J ;: .re^-eiyed ..signal y. For a. JDS -CDMA 
■ : -fX s ^ em -Vith repe^ output of the 

_ tjaird^ord^er ,Volterra_e^ ^ 

- ~ V-i*/v~r "~ ~ 

n-\ am 

. j. fl*0 cM) 

2..0 where _ the ^h cp-ef f icients are .estimated, usixig a. , linear 
n vt^^igue .such as MMSE , with ^.. the^_.autpcor^elat.ion and 
^ crosscorrelation .matrices, repl^ced -T w^th... the ^ autocorrelation 
; - .M^*. c^rpsscorrelatipn matrices , pf._ ,t^ ; j?Qwer ^se k rie,s... expanded 
matrices . "... 
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.-b.-:-..: ; j. -v. -; 

1. A direct sequence code division multiple access 
receiver comprising an adaptive filter controlled by an 

% '.adaptive- algorithm for r -filtering - . data which has been 
S "multiplied ,by -_a\: spreading r code . .and -filtered by a s channel 
filter, the adaptive filter having a length appropriate to 
-model the- -inverse... of., the channel filter, an d a .multiuser 
detector .operating.. on- the;.outputoOf the adaptive; filter . 

2. A receiver according; - to ;,cl t aim ■ 1 ; -wherein the 
10 algorithm is trained using the signal of a desired user. 

r. •; 3 : A.- receiver according to' daim l or-,2, wherein the 

algorithm, is ztra'ined. musing v-a composite ^signal from : more, than 
one user. 

; . : s*4 .if A ; receiver according to* claim. Tl.^.:2.-pr* 3 ,"- wherein 

15 the multiuser detector is of the minimum mean - squared" error 
type. 

5. A receiver according to claim 1, 2 or 3, wherein 
the multiuser detector is of the zero forcing 
(decorrelating) type. 

20 6. A receiver according to claim 1, 2 or 3 , wherein 

the multiuser detector is of the Volterra type. 

7. A receiver according to claim 1, 2 or 3 , wherein 
the multiuser detector is of the Radial Basis Function type. 

8. A receiver according to claim 1, 2 or 3 , wherein 
25 the multiuser detector is of the cancellation type. 

9. A receiver according to claim 1, 2 or 3 , wherein 
the multiuser detector is of the near optimum decoding type. 
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10 . A receiver According to any preceding claim, 
wherein the algorithm can comprises the least mean squares 
"algorithm. » ■ •_-.-> e.jc ~ " r . ; - - ..t- . fc 

- liv t A receiver- 'acco'rdfng -to any:one of . claims . 1_. to 9, 
5 wherein'fehe algorithm -comprises -She recursive: least squares 
algorithm. 1 " * - V- " r.-«r w ;. r/:: - - .:- ; i - 

' v - : : 1»2: — A- receiver" according: ft o^^n^r-^one vof .claims:, 1. to 9, 
wherein the algorithm comprises the fast a-posterion error 
sequential- technique k alg"6rithM^ r^v-.?:-. ; 

10 13 . A receiver according to any one of claims 1 to 9, 

wherein the" algorithm 1 a^ompri'se's " ^he^ " Stabilised fast a- 
J " * - posteriori e-3?rl5r --sequential technique":- 'algorithm.: \ sr^-g % 

14. A receiver according to claim 12 or 13, wherein 
*"■" ' said -algorithm ris . used in.:coinbifiautuioh:;w±th* the :Fast Newton 
15 -algorithm y^^rr r,jr,:j.\ - r • ;. z-zzzzj-i: r-a •/ : :r ;v: ? .: 
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